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Abstract 

Existence of discrepancies between experimental data and Monte Carlo predic¬ 
tions for angular distributions in four-jet production via e + e _ annihilation has 
been known for some time at LEP1. As such QCD processes constitute signifi¬ 
cant backgrounds to e + e _ —> W + W~ —> 4-jet production at LEP2, we consider 
the possibility that an erroneous modelling of the helicity structure of the final- 
state partons could affect the accuracy of experimental measurements of the W ± 
boson parameters. 
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1 Introduction and motivation 


The partonic reactions e + e~ —> qqgg and e + e~ —> qqQQ (see Fig. 1) constitute the dominant 
(lowest-order) component of the annihilation rate of electrons and positrons into four jets. 
They represent a benchmark process in QCD phenomenology at present and future e + e~ 
colliders for several reasons. 

First, they provide a test of the underlying SU ( Nc ) colour symmetry of the strong inter¬ 
actions between quarks and gluons, as the cross section is sensitive to the three fundamental 
colour factors C A , Cf (the Casimir operators of the fundamental and adjoint representations 
of the gauge group SU(Nc), respectively) and Tp (the normalisation of the generators of 
the fundamental representation). These are determined from the SU(Nq) generators ( T a )ij 
and the structure constants f abc via the relations 

J2f abc f abd * = 5 cd C A , £CrT°% = SijC F , Tr[T“T bt ] = 5 ab T F , (1) 

a,b & 

where a,b, ..(i, j ,...) represent gluon(quark) colour indices, yielding 

= N c , C F = (N 2 - 1)/2N C , T F = 1/2. (2) 

In QCD (i.e. Nc = 3) one has C A = 3 and C F = 4/3. The factors C A , C F and Tp represent 
the relative strength of the couplings (squared) of the processes q —> qg, g —> gg and g —> qq : 
respectively (see for example Ref. [1|]). 

Second, under the assumption that SU(N C = 3) is indeed the gauge group of the theory, 
a measurement of these colour factors (in particular of T fi = NpTp) can be converted into a 
constraint on the value of the number of quark flavours active at the energy scale at which 
the e + e _ annihilation takes place. For example, Np would be increased (by approximately 
3) from its (N F = 5) Standard Model (SM) value at LEP by an additional New Physics 
contribution from very light gluinos g [|2J produced in the process e + e _ —> qqgg via a g* —> gg 
splitting ||. 

Third, QCD four-jet production forms a significant background to e + e~ —* W + W~ —» 
4 jets at LEP2 []4[, for both the ‘threshold’ (total cross section) and the ‘direct recon¬ 
struction’ ( W ± line shape) methods employed in the measurement of M w , even after the 
implementation of typical W ± selection cuts [f|. Indeed, although biased by systematic 
uncertainties due to relatively unknown ‘Bose-Einstcin’ 0 and ‘colour interconnection’ 0 
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effects, the fully hadronic signature of the two W ± s still represents to date the experimen¬ 
tally preferred decay channel, because kinematic constraints can tighten the precision of the 
M\y measurement. 

Fourth, four-jet events represent a serious background in the search for Higgs particles 
at LEP2 [g and beyond (e.g. at the Next Linear Collider ||) JTT| , both in the standard 
electroweak (EW) theory and in possible extensions. In the SM, the dominant Higgs pro¬ 
duction channel proceeds via the ‘bremsstrahlung’ process e + e _ —> Z(j) JTT] , followed by the 
hadronic decays </> —*• bb. In this respect, it should be recalled that Z bosons have a 70% 
branching ratio into a pair of jets. Such arguments also apply to the case of the light h and 
heavy H Higgs scalars of the Minimal Supersymmetric Standard Model (MSSM) [|T2|. In 
addition, in the MSSM four-jet signatures arise naturally from the decays of the h scalar 
and the A pseudoscalar produced in pairs via e + e _ —*• hA 1 12; . 


Fifth, the four-jet channel could well be the footprint of new sparticles at LEP2 or be¬ 
yond, via the production of two resonances both decaying hadronically. In this context, 
it should be recalled that the production of several different types of supersymmetric par¬ 
ticles was advocated some years ago as an explanation of the apparent excess of four-jet 
events as recorded by the ALEPH Collaboration jnj. For example, pairs of sneutrinos [)I4| , 
squarks (in particular sbottoms |16]), charginos and neutralinos 0 , 0 ,m m and also 
selectrons [0] can yield four-jet signatures. Technicolour was also suggested as a possible 
explanation of the ALEPH events ^2 


For all of the above reasons it is clearly very important that four-jet events are correctly 
implemented in the Monte Carlo (MC) programs which are widely used in phenomenological 
studies of hadron production at e + e“ colliders (e.g. HERWIG 123)1 , JETSET/PYTHIA pdf . 
ARIADNE [0], etc.). In this connection, it is rather worrying that certain aspects of four-jet 
production are apparently not well described by the standard l O(a s ) ME + parton shower 
(PS)’ MC programs. For example, four-jet studies performed by the ALEPH Collaboration 


at LEP1 [0] have revealed a significant disagreement (up to 20% in some regions) between 
data and MCs for distributions in the four-jet angular variables: (i) the Bengtsson-Zerwas 
angle \bz ; (ii) the Korner-Schierholz-Willrodt angle <Lksw; (hi) the (modified) Nachtmann- 
Reiter angle (iv) the angle between the two least energetic jets # 34 . This is well 

illustrated by Fig. 2 of Ref. [|(| . 

The above angular quantities are defined in terms of the four-momenta of the particles 
in the final state and are particularly sensitive to the flavour (fermionic and bosonic) com- 
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position of four-jet events, in the sense that the different spins carried by quarks and gluons 
induce very different relative orientations between the two planes defined by the directions 
of the primary (i.e. from the 7*, Z'b-splitting) and secondary (i.e. from the virtual gluon 
splitting into quarks or gluons) parton pairsf]. For this reason, the angles Xbz, $ksw, $nr 
and 0 34 have profitably been used to fit the theoretical predictions to the data in order 
to measure the three colour factors Ca, Cf and T F . It can therefore be argued that the 
disagreement might be evidence that the MCs do not provide a correct description of the 
spin correlations among the various partons over the all phase space. Conversely, O(a^) 
ME programs (for example, the ‘(P(a^) ME + string fragmentation model’ as implemented 
in JETSET, i.e. with no parton shower) yield a much better angular description of four-jet 
final states (see Fig. 3 of ppj] and p^ ). Indeed, all spin correlations are naturally included 
in a full matrix element calculation (perturbative QCD predicts very specific orientations 
among the four hnal-state partons, see for example Ref. p|), but they are not necessarily 
present in a PS emulation of the four-jet final state, though their availability (in the infrared 
limit, where soft and collinear correlations can be factorised in analytic form) is a feature 
of some of the above mentioned MC codes (for example, for the HERWIG implementation, 
see Refs. |3T|). 

It might appear, therefore, that 0(a^) ME programs are the preferred tool for analysing 
QCD four-jet production at LEP (and beyond). However, the problem here is that such 
ME models contain ‘ad-hoc’ hadronisation which is adjusted to give a good description of 
some LEP1 data |26| (see also Ref. f32|) but cannot be reliably extrapolated to higher (e.g. 
LEP2) energies. Furthermore, it is a well known fact that their description of the sub-jet 
structure is very poor (see Ref. |h| for an overview). 

The deficiencies of these various approaches have been known for some time. Indeed, 
a general concern about, on the one hand, the inability of the standard MC programs to 
reproduce accurately four-jet angular quantities typical of parton level QCD and, on the 
other hand, the limitations of the ME implementation in quantifying hadronic multiplic¬ 
ities, already emerged at the time of the CERN LEP2 Workshop, see Ref. |I4|. In this 

3 To be more precise, the angle $34 was originally introduced to discriminate between the double 
bremsstrahlung (Fig. la b) and triple-gluon-vertex (Fig. lc) diagrams of the qqgg subprocess (see Ref. p7|), 
in order to provide direct proof of the non-Abelian structure of QCD. In fact, the variable $34 is more sen¬ 
sitive to the gluon propagator singularity in the graphs of Fig. lc, which is of course absent in those of 
Fig. la-b (see also the discussion in Sect. 3 below), than it is to spin correlations among partons. 
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respect, a key point not addressed there and that in our opinion deserves urgent attention 
is to assess whether or not also the angular behaviour of four-jet QCD events that survive 
W + W~ selection criteria is poorly modelled by the standard parton showers, as well as to 
quantify the sort of differences that one should expect from the two approaches. While it 
is the purpose of this study to specifically address this question, we also attempt here to 
understand the source of the disagreement between the PS and the ME predictions for the 
above angular variables. 

The plan of the paper is as follows. In the next Section we describe our method of 
computing four-jet rates using perturbative QCD. In Section 3 we present our results, and 
Section 4 contains a brief summary and some conclusions. 


e^e 


qqgg 


and 


e^e 


qqQQ events 


In order to generate the exact leading-order (LO) QCD predictions for the four-parton 
processes we use the FORTRAN codes already exploited in Refs. [|35], |36], |37|. Such programs 
include not only all the angular correlations between primary and secondary partons, but 
also all quark masses and both the 7 * and Z W intermediate contributions. As numerical 
inputs for the electroweak parameters we take a em = 1/128 and sin 2 9 W = 0.23, while for 
the Z boson mass and width we use Mz = 91.1 GeV and F z = 2.5 GeV, respectively. For 
the W ± boson, we use Mw = 80.23 GeV and F w ± = 2.08 GeV. The u, d and s quarks are 
taken to be massless, while for c and b we take m c = 1.35 GeV and nib = 4.95 GeV. The 
strong coupling constant a s has been computed at two-loop order (NLO), for N F = 5 active 
flavours, with Aq CD = 190 MeV and scale choice fi = >/s, where yfs is the centre-of-mass 
(CM) energy. 

In order to define a four-jet sample we need to introduce a jet algorithm. In the present 
study we will use the so-called ‘Durham’ (D) oncQ |39[], which is widely used at LEP and 
SLD, based on the ‘measure’ 


n 2 min (E?, Ej) (1 — cos 6G) 
Vij = --—:-• 


(3) 


In the above equation E t and Ej are the energies of the particles i and j, and dij is 
their relative angle (with i 7 ^ j = 1,.. .4). In our tree-level ME calculations, the four-jet 

4 See Ref. for a survey of jet clustering algorithms and their properties. 
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cross section for a given algorithm is simply equal to the four-parton cross section with a 
cut i/ij > y CVLt on all possible parton pairs (ij). In contrast, in our parton-shower studies, 
we need to define a jet clustering procedure as well. In these cases, we have adopted the 
so-called E-scheme, so that the four-momentum of a new cluster (or pseudoparticle) k is 
found from its constituents i and j by simple addition, p k = p t +Pj. The joining procedure 
is repeated until all pairs of clusters have a separation above y cut . This final set of clusters 
is called jets. 

At this point, we should further mention that we have also tried out three other jet 
clustering schemes, such as the ‘Jade’ |4(J, ‘Geneva’ |44|] and ‘Cambridge’ [f42] ones. Indeed, 
we have verified that none of our conclusions depends on the details of the jet algorithm 
used. 


For completeness, we recall the definitions of the four angular variables [43], [44], [45], [46 
introduced in the previous Section. One first orders the jets in energy such that 
Ei > E 2 > E 3 > Ei. This way, one can identify the two most energetic jets with those 
originated by the quarks produced in the Z decay, and the other two with those produced 
in the g* splitting (although such a procedure can be affected by systematic errors due to 
possible mis-assignments, see [f47] , f48f). Then, in terms of the three-momenta pi,...p 4 of 
the energy-ordered jets 1,... 4, the angles 6^ R , xbz and 634 are defined as 


6 nr = 4(pi ~P2,P3 -Pa), 
Xbz = Z(pl x p 2 j p 3 x p 4 ), 


(4) 

(5) 


and 


#34 = Z(p 3 ,p 4 ). ( 6 ) 

For <1>ksw we actually use the ‘modified’ definition proposed in Ref. |48|] (hereafter denoted 
by < h RSW ), which is more sensitive than the original to the flavour composition (see Refs. [44, 
PH). Thus, in events for which 

\pi+Ps\ > \pi+p4 (7) 

we define 

$ksw = Z(pi x p 3 ,p 2 x p 4 ), (8) 

whereas in the opposite case we dehne $ksw with P 3 and p 4 interchanged^. 

5 Note that the definition given in Eqs. (§)-(!) is equivalent to the original d> K sw angle Q in events 
where the thrust axis is directed along pj + p 3 or pi + P 4 . 
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3 Results 


Our results are presented in Figs. 2-8 and in Tab. I. They include no Initial State Radiation 
(ISR) at any of the energies, neither in the ME nor in the PS calculations. This is done 
in order to simplify the discussion, as we have verified that the relative behaviours of the 
two implementations are unaffected by the inclusion of such emission. In addition, notice 
that in the course of our analysis we will at times modify the value of the jet resolution 
parameter y cut adopted to define the four-jet sample. (In particular, the algorithm and the 
resolution used in Fig. 3 are the same as those of the ALEPH study |^ 6 |, enabling us to 
make quantitative comparisons of our findings with the results given there.) This will help 
to illustrate the generality of the phenomenology that we will describe and, in particular, 
will help to distinguish properties of the underlying matrix elements from artifacts of the 
jet definition. 

In the Introduction we argued that the disagreement between data and MCs could be due 
to a lack of adequate spin correlations among the partons generated in the phenomenological 
programs. In order to illustrate the sensitivity of the four angular variables to the hclicity 
structure of the 7 *, decay products we plot in Fig. 2 the distributions in xbz, $ksW’ $nr 
and 034 for the following subprocesses: (i) the four-quark component as obtained from the 
exact ME (tree-level) calculation; (ii) the same in a PS-like model obtained by averaging over 
the hclicities of the virtual gluon in Fig. Id (which removes the spin correlations between 
the two quark pairs in the final state); (iii) the triple-gluon-vertex component of two-quark- 
two-gluon events from the QCD ME (diagrams in Fig. lc); (iv) the same in a PS-like model 
again simulated by averaging over the virtual gluon spins. Note that we have not presented 
similar distributions for the case of the double bremsstrahlung component of the two-quark- 
two-gluon ME (see Fig. la-b), for two reasons^. First, it is strongly suppressed (by roughly 
a factor Ca/Cf = 9/4 |2iJ) with respect to the triplc-gluon-vertex contribution. Second, it 
is dominated by two infrared (i.e. soft and collinear) splittings of the virtual quarks into 
quark-gluon pairs, a dynamics which is well described by the PS implementation of the MC 
programs (as evidenced by their success in describing the phenomenology of e + e _ —> qqg 
production |29|]). In general, its angular behaviour is similar to that of the non-Abelian 
diagrams. 


6 See Ref. |4S] for some typical distributions in the above angles for all the components of the e + e 
parton process. 
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We emphasise that the results displayed in Fig. 2 (and also Fig. 4 below) should be 
regarded more as a useful exercise for understanding the underlying angular correlation 
properties of four-parton events in a ‘toy model’, rather than as a quantitatively reliable 
estimate of the relative behaviour of the exact QCD matrix elements vs. the MC parton 
shower description. In fact, two aspects should be noticed. First, the procedure of rewriting 
the diagrams in Fig. lc and d as the product of the ME for e + e _ —> qqg* times those of the 
decays g* —► gg and g* —> QQ, respectively, and of averaging over the spins of the intermedi¬ 
ate state (so that the secondary g* —> gg and g* —> QQ splittings are azimuthally symmetric 
about the virtual gluon direction), differs in several respects from that implemented in the 
MC programs. For example, here the large-angle splitting is correctly reproduced whereas 
in MC programs only the collinear part is correctly represented. In contrast, in our example, 
the soft splitting is described only at lowest order whereas in the parton cascade higher or¬ 
der (logarithmically) enhanced terms are also included. Furthermore, as already mentioned 
in the Introduction, some of the most sophisticated MC programs do include azimuthal 
and polar correlations in the parton branching (in the soft and collinear limit), so that our 
toy model of PS should be further regarded as an extreme condition, and for this reason 
particularly useful for our purposes. Secondly, our procedure of separating the Abelian 
and non-Abelian components in the two-quark-two-gluon partonic MEs by simply retaining 
the amplitudes associated with Fig. la-b and lc respectively, and neglecting their colour 
structure, is clearly non-gauge-invariant. However, we believe our distributions approximate 
quite accurately the true phenomenology of the four angular variables. In the first case, 
we do not expect the inclusion of large angle dynamics nor the neglect of logarithmically 
enhanced soft radiation to modify significantly the pattern of the angular correlations in 
the virtual gluon splitting. In the second respect, we stress that we have normalised all 
distributions to unity, since gauge invariance is more likely to affect the overall rates of the 
various components we have separated, rather than their angular shapes. We are therefore 
confident that the distributions shown in Fig. 2 (and 4) are an accurate representation of the 
behaviour produced by the corresponding diagrams in the full amplitude squared, and we 
can certainly use them as a guide to eventually pin-pointing the source of the discrepancies 
revealed in Ref. [pG[| . 

By looking at the Xbz distribution in Fig. 2, we see a clear peaking of the qqQQ com¬ 
ponent of the QCD ME around 90°, indicating the preference for the plane of the two 
secondary quark jets to be orthogonal to the plane of the two primary ones. In contrast, the 
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two secondary gluons from the triple-gluon vertex component prefer to be produced in the 
plane of the primary quark-antiquark pair (see also Fig. 9 of Ref. |I9|]). If in our toy model 
we switch off the angular correlations for the qqQQ final state we obtain a distribution that 
is significantly flatter, as expected. The ‘decorrelated’ version of the qqgg matrix element 
only including the triple-gluon-vertex diagrams gives a distribution similar to that of the 
decorrelated qqQQ model. The situation for the other three angular variables presented in 
Fig. 2, 'Fkswj ^nr and 0 34 , is very much the same as for xbz- That is, the two decorrelated 
versions of the QCD MEs are always significantly different from the ‘correct’ QCD ME 
prediction^. 

Therefore, Fig. 2 clearly confirms the usefulness of the four angles in testing the un¬ 
derlying structure of the QCD processes in the presence of the full ME spin correlations. 
Conversely, if the described spin dynamics is not taken into account exactly, then their 
discriminating power between qqQQ and qqgg events could be significantly reduced. This 
should be especially borne in mind when contemplating high-precision QCD analyses such as 
those aiming to determine the Ca, C f and T F colour factors (and/or the possible existence 
of light gluinos). 

For completeness, and in order to justify the confidence we have in our parton-levcl 
results, we also present in Fig. 3 the ratios between our ME and the PYTHIA results 
for the four angular distributions as obtained by adopting the same jet algorithm and 
resolution used in Ref. [25]. Here, in the MC implementation, partons are clustered before 
hadronisation until four jets are left. In addition, among all possible configurations of the 
latter, we only retain those for which the minimum of the measures is greater than y cut . 


The aim is to reproduce here the salient features of Fig. 2 of Ref. |s5|, our ME outputs 
mimicking the data of that experimental study. In this respect, note that in Fig. 2 of 
Ref. |25| a full simulation (DYMU + JETSET + ALEPH detector) was performed and 
the angular variables were computed at hadron level. However, given the good agreement 
between the two figures, we believe that the main conclusions that we will draw from the 
forthcoming studies are unaffected by our simplification. We will now set aside further 

7 Note that in Fig. 2 of Ref |2(| the distributions in the Bengtsson-Zerwas and (modified) Nachtmann- 
Reiter angles are ‘symmetrised’ by plotting with respect to |cos,\bz| and |cos0^ R |, respectively. This is 
because there exists an arbitrariness in the choice of the sign of the cosine of these angles. In contrast, in 
our Figs. 2 and 4 we have plotted the full angular range of the two variables, in order to gain more insight 
into the underlying problem of the parton shower implementation. Our distributions are not symmetric in 
themselves, since energy ordering of the jets is used (see for example the discussion in Ref. Q). 














considerations of the subject of precision QCD analyses and move on to another aspect of 
four-jet phenomenology. 

We consider here the more topical issue of W + W~ production and hadronic decays at 
LEP2. We have already mentioned that the QCD processes e + e _ —> 7 *, > qqgg, qqQQ 

constitute significant backgrounds to W + W~ —► 4-jet production at LEP2, up to 20% of 
the signal depending on the event selection criteria used. Thus, it is crucial that the MC 
generators used in the experimental studies are able to correctly describe the salient features 
(that is, the overall rate and the differential distributions) of the QCD four-jet final states 
that are compatible with the W + W~ —> 4-jet kinematics. 

A preliminary study of this problem was presented in Ref. |49|]. In particular, it was 
investigated there whether QCD events which pass the W + W~ event selection populate the 
ranges of the four-jet angular variables where the MCs have been shown to have problems 


describing the LEP1 data. Fig. 8 of Ref. |49|, obtained from simulations of QCD and W + W~ 
events at 172 GeV binned according to the four angular variables listed above, seemed to 
suggest that this is indeed the case. For example, using rather generic W + W~ selection 
cuts, it was shown that the W + W~ signal and the QCD background do cluster in the same 
angular regions. In particular, from Fig. 2 of Ref. IH and Fig. 8 of Ref. [f49| 
identify the following ranges of concern 


one can 


| cosxbzI, I cos^nrI < 0.7, | cos4>k SW |, I c °s 6> 34 1 < 0.5. 


(9) 


Using the QCD predictions taken from PYTHIA and the W + W~ predictions from KO- 
RALW |50|, one finds that at LEP2 approximately half the QCD events have, for example, 

| cos 6 * 34 1 < 0.5, a region largely populated by the signal and where the JETSET PS fails to 


describe the LEP1 QCD data by up to 15% |28 


It is therefore extremely important to investigate further such disagreement, particularly 
at LEP2 energies. For example, it is not impossible that W + W~ selection cuts might 
somehow affect the relative performances of the ME and PS implementations in describing 
the spin correlations. In particular, the constraints devised to disentangle the W + W~ signal 
at LEP2 could be such that the decorrelated MEs (and therefore also the MCs that they are 
supposed to emulate) approximate the exact results of the QCD MEs at LEP2 to a much 
better extent than they do at LEP1, thus reducing the severity of the problem. 

Unfortunately, this possibility does not appear to be true, as can be deduced from Fig. 4. 
This shows the same distributions as Fig. 2, but now at y / s = 172 GeV. In addition to the 
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processes studied at a/s = Mz we now include the corresponding angular distributions 
obtained at parton level from the e + e~ —> W + W~ —> qqQQ ME introduced and described 


in Ref. f5T]]. As a ‘typical’ W + W~ selection cut we have required M %] — Mw\ < 10 GeV 
for at least two pairs of partons (■ ij ). We see from Fig. 4 that, in the regions populated by 
W + W~ events, the differences in the shapes of the ME and PS-like implementations are still 
significant for xbz and <&k SW , somewhat reduced for 6 / R , and almost non-existent for $ 34 . 
However, before claiming that 6^4 has become a ‘safe’ quantity at LEP2, we should recall 
that, unlike the other variables, the angle between the two least energetic jets (generally, 
those produced in the virtual gluon splitting or by double gluon bremsstrahlung) is strongly 
dependent on the choice of the jet algorithm and/or the resolution parameter y cnt (i.e. on 
the higher-order terms of the perturbative expansion |)2| not included here). As we have 
already mentioned, it is a direct measure of the singularity which appears when the two 
secondary partons become collinear. (This singularity is of course absent for the W + W~ 
distribution, hence its flatness.) In other words, #34 is a rather ‘infrared unstable’ variable 
that, in our opinion, is not a reliable indicator of the underlying properties of the QCD 
four-jet matrix elements. 

Finally we turn to a comparison of the angular distributions at LEP2 obtained from the 
O(a^) QCD MEs and from two of the most widely used MC event generators, the up-to-date 
HERWIG 5.9 and PYTHIA 6.1, both of these using the l O{a s ) ME + PS’ approximation. 
We require them to reconstruct exactly four jets from the parton level before hadronisation, 
using the Durham jet finder with y cut = 0.002. Fig. 5a shows the distributions in the four 
angular variables at the LEP2 energy of 172 GeV with the usual W + W~ selection cut 
| Mij — M w | < 10 GeV enforced. 

Even assuming that the overall normalisation is the same in all cases (the distributions 
in Fig. 5a are all normalised to unity), we see clear discrepancies between the ME (solid 
histogram) and PS (dashed and dotted histograms) predictions, although apparently less 
pronounced than at LEP1 (see Refs. |26|, [28|]). Furthermore, the discrepancies are significant 
in the regions where the W + W~ signal peaks, especially for the Bengtsson-Zerwas and 
the (modified) Korner-Schierholz-Willrodt angles, well in line with our previous findings 
in the toy model. It is also evident that a simple rescaling of the distributions in the 
W + W~ populated regions does not resolve the differences. To quantify the size of the 
effects seen in Fig. 5a, we reproduce in Tab. I (Erst line) the fraction of events found in 
various angular intervals (typically, around the maximum of the W + W~ distributions) for 
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the ME predictions and those from HERWIG and PYTHIA. For reference, we also include 
the event fractions for the W + W~ signal. From Tab. I we see that the differences between 
the QCD ME and the PS predictions can be up to about 10-15%|. In addition, HERWIG 
and PYTHIA behave rather similarly. Thus, we conclude that one should expect significant 
discrepancies in the shapes of the angular distributions predicted by the ME and PS models 
when performing W + W~ —> 4-jet analyses. 

It is more difficult to draw any firm conclusions about possible differences in overall 
normalisation. The tree-level 0(a ME four-jet rate is rather unstable against the effects 
of NLO QCD corrections (see Ref. |]52||), with A'-factors which can be of (9(100%). How¬ 
ever, the shape of the angular distributions appears to be much more stable, implicitly 
demonstrating the relevance of the non-infrared dynamics in determining the behaviour of 
the four-jet angles. For example at LEP1, in the Durham algorithm and for y = 0.008 (the 
same set-up as in Ref. j^6[), the typical A^-factor of the angular spectra normalised to unity 


always lies in the range (1 ±0.025) for all angular variables, see Figs. 1-4 of Ref. There¬ 
fore, NLO effects alone cannot account for the discrepancies seen in Fig. 2 of Ref. nor 
do we expect them to be the reason for the differences noted in Fig. 5a here. For reference, 
we note that the JETSET LEP1 four-jet rates of (2(J and |28j are obtained by using an 
acceptance-rejection algorithm such that the three-jet rate of the generated events matches 
that given by the 0(a s ) ME. 

As we have mentioned the interplay between LO and NLO results in ME calculations, 
we would like to make one further comment. So far, we have compared the MEs at LO 
against the PS approximations of the MC programs, the latter defined by selecting exactly 
n = 4 jets with > y CVLt (with i % j = 1, ...4). Indeed, an alternative procedure which can 
be adopted to perform the comparison is to include also PS events with n > 4 jets (above 
the resolution) and eventually cluster these into exactly n — 4 jets. In fact, although the 
difference between the two approaches is next-to-leading and higher orders (recall that our 
ME are LO only), we have just mentioned that NLO effects do not alter the lowest order 
shapes of the MEs significantly (so should also be the case for NNLO, etc.). Furthermore, 
it is well known that in many cases summing over all possible radiation inclusively gives 
smaller ME/PS corrections than does vetoing resolvable radiation. Therefore, we present 
in Fig. 5b the usual angular distributions, now calculated from the two MCs in the new 


8 We have verified that the same size of effects are obtained using other jet algorithms and resolution 
parameters. 
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fashion. Indeed, the differences between the ME and the PS results are still sizable, as can 
be seen by looking at the second lines of Tab. I. 

Therefore, at this point it seems clear that a PS implementation based on 0(a^) QCD 
calculations for the hard scattering is needed in order to perform studies of four-jet events 
such as those outlined so far. In this respect, we note that an L 0(a^) ME + parton shower + 
cluster fragmentation’ option, based on both the correct second-order ME dynamics (from 
Refs. |)3], |56fl ) of the partons and supplemented with the appropriate showering of the 
latter, thus avoiding the shortcomings of the ‘(T(a^) ME + fragmentation’ (i.e., without 
PS) approach, will soon be publicly available in the new HERWIG version 6.1 ||54|p. The 
dashed curves in Fig. 3, showing the ratios of our MEs’ to the HERWIG 6.1 l O(a f) 

ME + PS’ implementation, at the parton level, well illustrates this improvement, allowing 
for residual differences between the two rates due to possible mis-assignments of the jet 
algorithm reconstructing the HERWIG parton shower backwards to the original four-parton 
state. The improvement is very much evident in the case of the angles xbz, &ksw an d $nr- 
The only exception occurs for 0 34 , where differences between the two approaches can still 
be relevant, especially when cos #34 —► 1. However, this is not surprising, as we have already 
highlighted the sensitivity of this variable to higher-order effects in the collinear limit, which 
are indeed embodied in the PS evolution but not in the LO MEs (nor in the l O(ag) ME 
+ hadronisation’ model). Not surprisingly, the full NLO corrections to the shape of this 
angular distribution can be as large as 15%, again in the region 0056*34 —> 1, for example 
in the Cambridge scheme f53|. The appropriate treatment of such Sudakov effects due to 


soft-gluon emission |29[ are beyond the scope of the present study and will be addressed 
elsewhere [[54 . 


One might now wonder whether the above parton-level behaviours can survive the hadro¬ 
nisation process. We address this in Fig. 6 , where the usual four angular distributions are 
plotted, but now at the hadron level (including both charged and neutral tracks in the jet 
reconstruction and allowing for the decays of the heavy hadrons), as obtained by HERWIG 
6.1, using both 0(a s ) and O(a^) initiated PSs. From Fig. 6 , it is evident that the soft dy¬ 
namics of the hadronisation stage does not remove the partonic differences between the two 
approaches. Furthermore, by comparing Fig. 5a to Fig. 6 , one can see that the discrepancies 
are of the same sort at the two QCD stages, being possibly larger at the hadron level in 


9 Some progress in the same direction as in HERWIG is being made also in the context of the PYTHIA 
environment [B3|. 


12 












the case of the Bengtsson-Zerwas angle. Therefore, it is more than plausible that ‘partonic’ 
angular effects are the source of the differences between the current PS MCs and the LEP1 
data. 

A very relevant issue which can be addressed at this point is whether such effects can 
influence the determination of physical parameters whose measurement depends on the 
estimation of the size and topology of the four-jet background from QCD. One important 
and topical example concerns the measurement of the W mass at LEP2, reconstructed from 
fully hadronic decays of W ± pairs, as discussed in the Introduction. However, in order 
to give quantitatively reliable estimates in this case, one should take into account various 
other non-trivial effects such as initial-state QED radiation, neutrinos that escape without 
detection, cracks in the detector acceptance, experimental resolution and efficiency, etc. In 
addition, we note that the various selection procedures for W + W~ —> 4-jet candidate events 
differ significantly from experiment to experiment. Therefore it is clear that a full study of 
the effect on the W mass can only be carried out within the context of a complete detector 
simulation, which is beyond the scope of this paper. 

Nevertheless, we can address this issue here at a rather more modest level, to gauge the 
size of the effect. Thus, we run HERWIG once more to produce hadronic final states via the 
signal process e + e _ —> W + W~ —> qqQQ (IPR0C=200, as in 5.9) and the background e + e _ —> 
7 *, Z* —> qqgg, qqQQ , the latter generated both via 0(a s ) (IPR0C=100, as in 5.9) and 0(a 2 s ) 
(IPR0C=600, new to 6.1) MEs. Again, we require exactly four jets to be reconstructed by 
means of the Durham scheme with y cut > 0.002 and such that M t j — My/\ < 10 GeV for at 
least two pairs of jets (ij ), but this time only within the angular regions typically populated 
by the signal. Following Fig. 5a, we have chosen the following constraints^]: 

l cos X-Bz| > 0.5, cos^ksw < — 0.5, |cos#n R | > 0.5 |cos# 34 | < 0.8. (10) 

Next, we plot an ‘average’ W ± mass, hereafter denoted by M ave , for each event generated. 
There are several reasons for choosing the average rather than the two individual masses. 
For example, the mis-assignment of one particle from one W ± to the other reduces the first 
mass and increases the second, leaving the average less affected than each separately. In 
fact in each event there are three possible jet pairings, each giving one potential average W 
mass. Of these three, we exclude the one where the two most energetic jets are paired with 

10 Given the remark in Footnote 7, recall that we ought to symmetrize the cuts in the Bengtsson-Zerwas 
and Nachtmann-Reiter angles. 
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each other, since kinematically this is seldom the correct combination. The plot in Fig. 7 
shows the shape of such a distribution for the three cases IPR0C=200 (solid line), IPR0C=100 
(dashed line) and IPR0C=600 (dotted line). From this figure the different shape of the two 
descriptions of the backgrounds is rather clear, particularly in the vicinity of the W ± mass 
resonance. 

As a further step in our study we bin (see upper part of Fig. 8 ) the average W ± mass 
spectra as obtained by summing the signal and background rates, each with the appropriate 
normalisation as computed by HERWIG, in the two possible cases, depending on whether 
the (D(a s ) (dashed line) or the 0(a 2 s ) (dotted line) MEs are used in the MC event generator 
to simulate the QCD noise. The ratio in each bin between the latter two distributions is 
shown in the lower plot of Fig. 8 (solid line). Indeed, the total spectra also differ significantly 
and the ratio between the two is not constant around M\y. 

To investigate the possible consequences of the results shown in Figs. 7-8 on the deter¬ 
mination of the W ± mass, we perform a MINUIT [|57j fit on the two total distributions, with 
a fitting function of the form 


f(m) = ci 


,2 „2 


CoC 


2^3 


(m 2 - cl ) 2 + c|c 2 


+ g{m) 


( 11 ) 


where the term g(m) is meant to simulate a smooth background due to mis-assigned jets 
induced by the clustering algorithm. For the latter, we adopt three different possible choices 


g(m) = 


0, 


c 4 + c 5 (m - c 2 ) + c 6 (m - c 2 ) 2 , 

1 +exp ((m—C5) /c 6 ) ’ 


( 12 ) 


that is, a null, a three-term polynomial and a smeared step function (motivated by the 
kinematical-limit shoulder at large masses). Note that in eq. (|TTD we have assumed a Breit- 
Wigner shape characterised by a peak height ci, a position c 2 and a width C 3 , corresponding 
to the normalisation, Mw and Ty y, respectively, of the two distributions in Fig. 8 . 

To first approximation, then, the difference between the two values of the coefficient c 2 
as obtained from fitting the two curves is a measure of the typical size of the systematic 
error that could be introduced in the experimental measurement of the W boson mass by a 
mis-modeling of the QCD four-jet background. We find this value to be around 10 MeV if 
we neglect altogether the intrinsic mis-assignment background or if we assume for the latter 
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a step function, whereas the difference somewhat decreases for a polynomial background, 
down to 2 MeV. Though not dramatically large, such numbers are nonetheless comparable 
to the expected size of the final systematic error on Mw from all other possible sources (as 
estimated for example at the LEP2 Workshop ||) of about 50 MeV. 

4 Summary and conclusions 

In this paper we have investigated discrepancies in the predictions for distributions of typical 
four-jet angular variables obtained from matrix element and parton shower calculations. 
Comparison of the latter with LEP1 data from ALEPH revealed differences of up to 20%, 
while the ME calculations are in good agreement with the same sample. By constructing 
toy parton-shower models, using matrix elements in which the spin correlations between 
the secondary (i.e. from the g* splitting) QQ and gg pairs and the primary (i.e. from the 
7 *, Z W decay) qq ones are switched off, we concluded that the source of the disagreement 
resides in the limited capability of the MC parton shower of reproducing the exact partonic 
spin correlations away from the infrared limit, if the ‘0(a s ) ME’ is used to generate the 
hard QCD scattering. 

However, the main concern of our analysis was the impact that similar effects could 
have on four-jet studies at LEP2, in the context of W + W~ —> 4-jet phenomenology and 
W ± mass determinations. The preliminary exercises that we have carried out in this paper 
show that the shape of the angular distributions is not reproduced accurately by the MC 
parton showers at LEP2 either. Furthermore, these discrepancies survive even after the 
W + W~ selection criteria are implemented as well as the hadronisation phenomenon. 

It is possible, therefore, that at present the contribution of the QCD four-jet background 
is not being correctly simulated in candidate W + W~ samples, and this could constitute a 
non-negligible source of error in the M\y determination at LEP2. We regard this as the 
most important conclusion of our paper, based on the fact that, after having performed a 
simple exercise to assess the size of such systematic error, we have found it to be as large as 
10 MeV. This is of the same order as the foreseen precision on M\y at the end of the LEP2 
runs. However, we stress that our estimate should not be taken too literally, since our study 
lacks many of the necessary ingredients of a complete simulation, such as detector effects, 
ISR, a more detailed background analysis (including all four-quark EW channels), etc., all 
of which is well beyond the scope of this paper. 
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Nonetheless, it is evident that it is of crucial importance that a dedicated implementation 
based on O(a^) QCD calculations is used in the MC simulations of four-jet samples of the 
type described in this paper. MC programs are now being upgraded to incorporate such 
second-order dynamics and we urge the experimental collaborations to adopt them in their 
analyses. 
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Table Caption 


[I] Fraction of four-jet events in various angular intervals around the maxima of the 
W + W~ distributions, as predicted by the: (i) MEs’ (qqgg + qqQQ) , (ii) l O(a s ) 

ME + PS’ (HERWIG, IPR0C=100), (hi) l O(a s ) ME + PS’ (PYTH1A, ISUB=1); and 
also, for reference, the l O(a®) ME’ for e + e _ —> W + W~ —> qqQQ. All rates are 
obtained at parton level. The following cut has been implemented: | My — Mw \ < 10 
GeV on at least two pairs of partons (ij). The CM energy is y/s = 172 GeV. Results 
are shown for the Durham algorithm, with y cut = 0.002. The first line refers to the 
case in which exactly n — 4 jets are reconstructed from the parton shower, while the 
second corresponds to n > 4, with the additional jets in the latter eventually clustered 
into four. (All n final tracks are above cut-off in both cases.) 



Figure Captions 


[1] Feynman diagrams contributing in lowest order to QCD four-jet production in e + e“ 
annihilation: (a,b) double-gluon-bremsstrahlung, (c) triple-gluon-vertex and (d) four- 
quark subprocesses. 

[2] Differential distributions for the following variables: (top-left) xbz, (top-right) <h RSW , 
(bottom-left) cos d^ R and (bottom-right) cos #34 for the following subprocesses: four- 
quark in QCD (solid), four-quark in a PS-like model (dashed), triple-gluon in QCD 
(dotted), triple-gluon in a PS-like model (dot-dashed). The CM energy is yfs = Mz- 
All distributions are normalised to unity. Results are shown for the Durham algorithm, 
with y cut = 0 . 002 . 

[3] The ratio between the ME and PYTHIA implementation (solid) of the differential 
distributions for the following variables: (top-left) xbz, (top-right) f^w, (bottom- 
left) cos 6 ^ and (bottom-right) cos # 34 . The default spin correlation matrices of the 
parton cascade in the MC program have been switched on. All distributions are at 
parton level. The CM energy is y/s = Mz- All distributions are normalised to unity. 
Results are shown for the Durham algorithm, with y cllt = 0.008. Exactly four jets are 
required to be reconstructed from the PYTHIA cascade. The corresponding ratio in 
case of HERWIG 6.1 is also shown (dashed): see later on. 

[4] Differential distributions for the following variables: (top-left) xbz, (top-right) <I> RSW , 
(bottom-left) cos and (bottom-right) cos #34 for the following subprocesses: four- 
quark in QCD (solid), four-quark in a PS-like model (dashed), triple-gluon in QCD 
(dotted), triple-gluon in a PS-like model (dot-dashed) and W + W~ (shaded). The CM 
energy is yfs = 172 GeV. All distributions are normalised to unity. The following cut 
has been implemented: | M tJ — Mw \ < 10 GeV on at least two pairs of partons (ij ). 
Results are shown for the Durham algorithm, with y cut = 0.002. 

[5] Differential distributions for the following variables: (top-left) xbz, (top-right) ^ksw; 
(bottom-left) cos^r and (bottom-right) cos 6*34 in case of l O(al) MEs’ (qqgg + qqQQ, 
solid histogram), l O(a s ) ME + PS’ (HERWIG, dashed histogram) and l O(a s ) ME + 
PS’ (PYTHIA, dotted histogram). The default spin correlation matrices of the parton 
cascade in the MC programs have been switched on. For reference, we also have 



included the same spectra as obtained from the ‘ 0(a °) ME’ for e + e _ —» W + W~ —> 
qqQQ (asterisk symbols). All distributions are at parton level. The following cut 
has been implemented: — M\y\ < 10 GeV on at least two pairs of partons (ij). 

The CM energy is y/s = 172 GeV. Normalisation is to unity. Results are shown for 
the Durham algorithm, with y CVLt = 0.002. (a) Exactly four jets are required to be 
reconstructed from the HERWIG and PYTHIA cascades, (b) At least four jets are 
required to be reconstructed from the HERWIG and PYTHIA cascades, eventually 
forced into exactly four. 


[ 6 ] Differential distributions for the following variables: (top-left) xbz, (top-right) ‘hxsw; 
(bottom-left) cos 6 ^ and (bottom-right) cos #34 in case of l O(a 2 ) ME + PS + Hadro- 
nisation’ (solid histogram) and c O(a s ) ME + PS + Hadronisation’ (dashed histogram), 
as obtained by using HERWIG. The default spin correlation matrices of the parton 
cascade in the two versions of the MC program have been switched on. The following 
cut has been implemented: M t] — Mw\ < 10 GeV on at least two pairs of jets (ij), 
as defined at the hadron level. The CM energy is yfs = 172 GeV. Normalisation is to 
unity. Results are shown for the Durham algorithm, with y cut = 0.002. Exactly four 
jets are required to be reconstructed. 


[7] Differential distributions in the average H /± -mass (as defined in the text) for the 
W + W~ signal (solid curve), the 0(a s ) (dashed curve) and the 0(a 2 s ) (dotted curve) 
backgrounds (from HERWIG). Other than the constraint — Mw\ < 10 GeV on at 


least two pairs of jets (ij), as defined at the hadron level, the angular cuts of eq. ( p~ 0 |) 
have also been implemented. The CM energy is y/s = 172 GeV. Normalisation is to 
unity. Results are shown for the Durham algorithm, with y cnt = 0.002. Exactly four 
jets are required to be reconstructed. 

[ 8 ] Differential distributions (top) in the average IFVmass (as defined in the text) for 
the sum of the W + W~ signal and the 0(a s )(0(a 2 )) background (from HERWIG): 
dashed(dotted) curve. Other than the constraint \M tJ — Mw \ < 10 GeV on at least 
two pairs of jets (ij), as defined at the hadron level, the angular cuts of eq. (|T0|) have 
also been implemented. The CM energy is yfs = 172 GeV. Normalisation is to the 
total rates. Results are shown for the Durham algorithm, with y cut = 0.002. Exactly 
four jets are required to be reconstructed. The bottom plot shows the ratio between 
the above two spectra. 
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